In-flight dissipation as a mechanism to suppress Fermi acceleration. 



O 

(N 

> 
O 



Q 

U 
d 



> 

OS 

o 



X 



Diego F. M. Oliveira and Marko Robnik 
CAMTP - Center For Applied Mathematics and Theoretical Physics 
University of Manbor - Krekova 2 - SI-2000 - Maribor - Slovenia. 
(Dated: November 2, 2010) 

Some dynamical properties of time-dependent driven elliptical-shaped billiard are studied. It was 
shown that for the conservative time-dependent dynamics the model exhibits the Fermi acceleration 
[Phys. Rev. Lett. 100,014103 (2008)]. On the other hand, it was observed that damping coefficients 
upon collisions suppress such phenomenon [Phys. Rev. Lett. 104, 224101 (2010)]. Here, we consider 
a dissipative model under the presence of in-flight dissipation due to a drag force which is assumed 
to be proportional to the square of the particle's velocity. Our results reinforce that dissipation 
leads to a phase transition from unlimited to limited energy growth. The behaviour of the average 
velocity is described using scaling arguments. 

PACS numbers: 05.45.Ac, 05.45.Pq 



Dissipative systems have attracted much attention 
during the last years since they can be used in order to 
explain different physical phenomena in different fields of 
science including atomic and molecular physics 0, Q , tur- 
bulent and fluid dynamics [3|-l5| , optics [6j, [D] , nanotech- 
nology d, IH , quantum and relativistic systems [l(| EJ ■ 
Different procedures can be used to describe such sys- 
tems. The billiard models are often considered since they 
can be easily described mathematically and can be real- 
ized experimentally in many different ways, for example, 
the microwave resonators initiated by H.-J. Stockmann in 
1990 [l2| and also superconducting microwave resonators 
[l3j|, quantum dots [HI], ultracold atoms [lj| and many 
others. From the mathematical point of view, a billiard 
is defined by a connected region Q C R D , with boundary 
dQ C R 0-1 which separates Q from its complement. If 
dQ = dQ(t) the system has a time-dependent boundary 
and it can exchange energy with the particle upon col- 
lision. In such a case, it is possible to investigate the 
phenomenon called Fermi acceleration, i.e., the unlim- 
ited energy growth [n| . According to Loskutov-Ryabov- 
Akinshin (LRA) conjecture [l7|, a chaotic component in 
the phase space with static boundary is a sufficient condi- 
tion to observe Fermi acceleration when a time dependent 
perturbation is introduced. Results that corroborate the 
validity of this conjecture include the time dependent oval 
billiard [ill EH > stadium billiard [2(| , Lorentz Gas [2l[ . 
Recently, it was shown even that a specific perturbation 
in the boundary of an elliptical billiard (integrable) leads 
to the unlimited energy growth 22[ . The separatix gives 



place to a chaotic layer and the particles can now ex- 
perience unlimited energy growth while diffusing in the 
chaotic layer. 

In this Letter, we will consider a dissipative ellipti- 
cal billiard with a periodically moving boundary which 
has been studied in the pioneering paper in 1996 [24j . 
Firstly, wc assume that the particles of mass m are im- 
mersed in a fluid. The dissipative drag force is consid- 
ered to be proportional to the square of the velocity of 
the particle, V . To obtain the equation that describes 
the velocity of the particle along its trajectory, we need 



to solve Newton's equation where m 
with the initial velocity \? n > and if is the coeffi- 
cient of the drag force. After we introduce the vari- 
ables rj'/m = rj, we obtain the velocity of the particle 

V„ 



as function of time as V p (t) = 



We de- 



1+| V n \n(t-t n ) 

scribed the model using a four dimensional and non lin- 
ear map T[9 n , a n , \ v n \,t n ) = (6 n +i, a n -|_i, 

; tn+l) 

where the dynamical variables are, respectively, the angu- 
lar position of the particle; the angle that the trajectory 
of the particle forms with the tangent line at the position 
of the collision; the absolute velocity of the particle; and 
the instant of the hit with the boundary. Figure [1] illus- 
trates the geometry of five successive collisions of the par- 
ticle with the time-dependent boundary. To obtain the 
map, we start with an initial condition (9 n , o; n , v n 
The Cartesian components of the boundary at the angu- 
lar position (9 n ,t n ) are 



X(6 n ,t n ) = [^ + C , sin(i„)]cos(6» r , 
Y(6 n ,t n ) = [B + Csm(t n )]sm(9 n 



(1) 
(2) 



where Aq and Bq are constants, thus, at any time t n we 
have elliptical shape. The control parameter < C < 
min(Ao, Bq) controls the amplitude of oscillation and 
9 € [0, 2ir) is a counterclockwise polar angle measured 
with respect to the positive horizontal axis. The an- 
gle between the tangent of the boundary at the position 
(X(9 n ), Y(9 n )) measured with respect to the horizon- 

where the expressions 

and Y'(9 ni t n ) = 



= arctan 



Y'(8 n ,t n ) 



tell lillG is i^n — cLLK>ban x.' (9 t ) 

for both X'(9 n ,t n ) = dX(9 n ,t n )/d9 n a ^ * y u n ,, 
dY(9 ni t n )/d9 n . Since the expressions for <f> n and a„ are 
known, the angle of the trajectory of the particle mea- 
sured with respect to the positive X-axis is (</>„ + a n ). 
Such information allows us to write the particle's veloc- 
ity vector as V n = \V^,(t)\[cos((f> n +a n )i + sm(<f> n +a n )j}. 
Where i and j denote the unity vectors with respect to 
the X and Y axis, respectively. The particle travels on a 
straight line until it hits the time dependent boundary. 
The position of the particle, as a function of time, for 
t > t n , is X p (t) = X(6 n , t n ) + r(t) cos(0„ +«„), Y p (t) = 



FIG. 1: Illustration of five collision with the time dependent 
boundary. The corresponding angles that describe the dynam- 
ics are also illustrated. 
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FIG. 2: Phase space for the static elliptical billiard. The 
control parameters: Ao — 2 and Bo = 1. 



Y(9 n , t n ) + r(t) sin(</>„ + a n ). Where the sub-index p de- 
notes that such coordinates correspond to the particle 
and r{i) = r\~ x ln[l + V n n(t — t n )], which is the displace- 
ment of the particle obtained from direct integration of 
dr(t)/dt = v p (t). The distance of the particle measured 
with respect to the origin of the coordinate system is 

given by R p {t) = Jx$(t)+Y*(t) and 9 P at X p (t), Y p (t) 



is 6 P — aTctan[Y p (t) / X p (t)]. Therefore, the angular 
position at (n + l) th collision of the particle with the 
boundary, i.e. n +i, is numerically obtained by solving 
the following equation R p {t) = ^/X 2 (6 p ,t) +Y 2 (6 p ,t). 
The time at (n + l) th collision is obtained evaluating 
t n+ \ = t = t n + t c , where t c is the time during the flight. 
To obtain the new velocity we should note that the ref- 
erential frame of the boundary is moving. Therefore, at 
the instant of collision, the following conditions must be 
obeyed 

V n +i • T n _|_i = T^„ • T n+ i , (3) 

Vn+l ■ ?$n+l = -fn • ^ n+l + if b (t n+1 ) ■ n+1 (,4) 

where the 1* and 7$ are the unitary tangent and normal 
vectors, respectively, and the velocity of the boundary 
Vb{t n +i) = C cos(t n+ i)[[cos(6 n +i)i + sm(9 n+1 )j}. Then 
we have 



^n+l\ = J(V n+ l ■ fn+l) 2 + ■ ^ n+l) 2 . (5) 



Finally, the angle cv n +i is written as 
a n +i — arctan 



V n+l - ~^n+l 
~^ n+l ■ ^ n+l 



(6) 



With this four dimensional mapping, we can explore 
the dynamics of the model. However, before considering 
the time-dependent model let us illustrate the behaviour 
of the phase space for the static boundary. Indeed, it 
is well known that the ellipse is an integrable billiard 
system, the product of the two angular momenta with 
respect to the foci being the integral of motion [2^, [26[ . 
Figured! shows the phase space for A = 2 and B = 1. 
We can see a large double island limited by a separatrix 
and a set of invariant spanning curves. We can observe 



two different kinds of behaviour separated by a separatrix 
(red curve), namely, rotators and librators. Librators 
consist of trajectories that are confined between the two 
foci and in the phase space are confined by the separatrix 
curve. On the other hand, rotators are trajectories near 
to the boundary exploring all the values of 8. In the 
phase space they are outside of the separatrix curve. 

As a part of our numerical results we shall discuss 
mainly the behaviour of the average velocity of the par- 
ticle. Two different procedures were applied in order to 
obtain the average velocity. Firstly, we evaluate the av- 
erage velocity over the orbit for a single initial condition 
and then over an ensemble of initial conditions. Hence, 
the average velocity is written as 



V = 



M n 

M ^ n + l ^ ■ 



(7) 



where the index i corresponds to a sample of an ensem- 
ble of initial conditions, M denotes the number of dif- 
ferent initial conditions. We have considered M = 200 
in our simulations. It was shown by Lenz et al. 0, [23| 
that when a driving perturbation is introduced into the 
system, opposite to the expectations, it presents a phe- 
nomenon known as Fermi acceleration or unlimited en- 
ergy gain. Such a behaviour happens because when the 
driving amplitude C ^ the separatix is replaced by a 
chaotic layer. A particle which starts its dynamics in a 
rotator orbit can change its dynamics to a librator and 
vice versa. The chaotic diffusion within the chaotic layer 
leads to unlimited energy growth. Figure [3] shows the be- 
havior of the average velocity as a function of the number 
of collisions. We have considered the conservative case 
where the drag coefficient is n — 0. For such a case the 
particle's velocity is V p (t) — v n andr(i) = |K n |(t— t n ). 
As one can see, all curves of the V behave quite similarly 
in the sense that: (a) for short n, the average velocity 
remains constant for a while, but eventually, (b) after 
a crossover, all the curves start growing with the same 
exponent. This is at variance with the result obtained 
by Lenz et al [22j ], perhaps because they were not con- 
sidering the average velocity for large enough values of 
n. 
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FIG. 3: Behaviour of Vvs.n for different initial velocities. 
The control parameters used were Aq = 2 , Bq = 1. We have 
" '•idered the coefficient of the drag force as r\ = 0. 
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We discuss now the effect of dissipation introduced via 
frictional force. To obtain the average velocity, we ran- 
domly choose t G [0,27r], G [0,2tt] and a G [0,tt]. We 
also fix the value of n = 10 ~ 3 . Additionally, in order to 
avoid the initial plateau we also have fixed initial veloc- 
ity as Vb = 10~ 5 . The model of collisional dissipation by 
Leonel and Bunimovich 29] is different from our model 
of in-flight dissipation due to the drag force in detail, but 
should behave similarly in the statistical sense (on the av- 
erage), especially in the chaotic regime, because then we 
have < V n +\ >=< V n > e~ vr , thus the effective damp- 
ing coefficient 5 = e~ ri<r> , where < r > is the mean free 
path of the particle. 

In Fig. @] (a) we show the behaviour of the average 
velocity as a function of the number of collisions for dif- 
ferent values of C. Note that for different values of C and 
for short n, the average velocity starts to grow and then 
it bends towards a regime of saturation for long enough 
values of n. It must be emphasized that different values 
of the parameter C generate different behaviors for short 
n. However, applying the transformation n — > nC 2 coa- 
lesces all the curves at short n, as shown in Fig. 0|b). For 
such a behaviour, we can also propose the following scal- 
ing hypotheses: (i) When n <C n x the average velocity 



V(nC\C) oc (nC A ) 



(8) 



where the exponent /3 is called the acceleration exponent, 
(ii) When n^> n x , the average velocity is described as 



V sat oc CP. 



(9) 



where 7 is the saturation exponent, (iii) The crossover 
from growth to the saturation is written as 



n x oc C z . 



(10) 



where z is called crossover exponent. 

These scaling hypotheses allow us to describe the av- 
erage velocity in terms of a scaling function of the type 




FIG. 4: (a) Behaviour of the average velocity as function of n 
for different values of the control parameter C. (b) Their ini- 
tial collapse after the transformation nC 2 . (c) Their collapse 
onto a single universal plot. 



where p and q are scaling exponents and A is a scaling 
factor. Moreover, p and q must be related to j3, 7 and 
z. Because A is a scaling factor, we can specify that 
\PnC 2 = 1, yielding 



V[nC 2 ,C] = (nC 2 ) 1/p V l [n~ q '- p C] , 



(12) 



where V 1 [{nC 2 )- q / p C] = V[l, (nC 2 )- q / p C] is assumed 
to be constant for n <C n x . Comparing Eqs. ([8j and 
(fl"2"l) . we obtain (3 — —1/p. Choosing now X q C — 1, we 
find that A = C" 1 / 9 and Eq. dTTJ) is given by 



V[nC 2 ,C] = C' 1/q V 2 [C~ p,q nC 2 



(13) 



V[nC 2 ,C] =XV[X p nC 2 1 \ q C] , 



(11) 



where V 2 [C- p/q nC 2 ] = V[C~ p / q nC 2 , 1] is assumed to be 
constant for n 3> n x . Comparing Eqs. (|9|) and (p~3|) . we 
obtain 7 = — 1/q [see Fig. [5] (a)]. A power law fitting 
in Fig. [5] gives us that [3 = 0.501(2). Such value was 
obtained from the range of C G [10 _3 ,2 x 10 -1 ]. Given 
the two values of the scaling factor A, one can easily con- 
clude that z = -g — 2 = —0.97(1), which is in excellent 
agreement with the value obtained numerically, as shown 
in Fig. [5] (b) . A confirmation of the initial hypotheses 
is made by a collapse of all the curves of Vvs.n onto a 
single and universal plot, as shown in Fig. |4] (c), show- 
ing that the system is scaling invariant under specific 
transformation. With this good collapse of all the curves 
of the average velocity and considering that the critical 
exponents are j3 = 0.5, 7 = 0.5 and z = —1, we can 
conclude that the time dependent driven elliptical bil- 
liard belongs to the same class of universality of the one 
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FIG. 5: (a) Plot of Vsat as function of the control parameter 
C. (b) Behaviour of the crossover number n x against C. 



dimensional Fermi-Ulam model [27| and the periodically 
corrugated waveguide [Hj]. The scaling can also be de- 
scribed in terms of is the coefficient of the drag force 77. 
In such a case V sat oc 7y-°- 521 (3) anc } Hx K ^-1.070(6) < ^ e 

have fixed C = 0.1 and Vq = 10~ 5 . Therefore, 77 — ^ im- 
plies that V sa t and n x both diverge, thus recovering the 



results for the conservative case, i.e., exhibiting Fermi 
acceleration. Additionally, our results reinforce that in- 
flight dissipation is a sufficient condition to suppress the 
phenomenon of Fermi acceleration like in the case of col- 
lisional dissipation 29] . 

As the concluding remark, we have studied some dy- 
namical properties of a time dependent driven elliptical 
billiard. In the static case the phase space is integrable 
where two kinds of trajectories are observed: rotator and 
librator. We have introduced time dependent perturba- 
tions and in-flight dissipation. We have observed that 
average velocity grows for small number of collision and 
then, after a crossover, it reaches a regime of saturation 
for large n. Thus we do not observe the unlimited en- 
ergy growth (Fermi acceleration). We have also studied 
the behaviour of the average velocity using scaling argu- 
ments. We have shown that there is a relation between 
the critical exponents 7, ft and z. Our scaling hypotheses 
are confirmed by a perfect collapse of all the curves onto 
a single universal plot. Additionally, we confirm that 
the two dimensional elliptical model belongs to the same 
class of universality of the Fermi-Ulam model (1-D) and 
the corrugated wave guide (1-D), for the range of control 
parameters studied. 
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